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_^ . Abstract 

In this paper we show how one can obtain simultaneous rational approximants for 
\^ , C(j(l) s-nd C(j(2) with a common denominator by means of Hermite-Pade approxima- 

r \ I tion using multiple little g-Jacobi polynomials and we show that properties of these 

• ■ rational approximants prove that 1, Cg(l)) Cg(2) are linearly independent over Q. In 

"^ , particular this implies that C(?(l) aiid C,q{2) are irrational. Furthermore we give an 

■ ' upper bound for the measure of irrationality. 
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^ ' 1 Introduction 

o 

'^ I There are a lot of things known about the irrationality of the Riemann zeta-function 

OO 1 

^ : at integer values s e {2,3,4,...}. It is known that C(2n) = (-l)"-^22"-i52„7r2"/(2n)!, 

where -82™ are Bernoulli numbers (which are rational), so it follows that C{2n) is an 
irrational number for n G Nq = {1,2,3,...}. It is also known that ({3) is irrational 
(Apery [2]) and that C{2n + 1) is irrational for an infinite number of n G Nq [3] [6] [21]. 
Furthermore at least one of the numbers C(5), C(7), C(9)) C(ll) is irrational [16] [20]. 

A possible way for a g-extension of the Riemann zeta-function is 

Us) = T.-r^^ . = 1,2,..., (1.1) 

n.=l -*- ^ 
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with g G C and |g| < 1. Then the hmit relations 

hm (1 - q)Xg{s) = {s-iy. C(s), 5 = 2, 3, ... , (1.2) 

kl<i 

hold [23]. Earlier it was shown that Cq{l) is irrational whenever q = 1/p, with p an integer 
greater than 1 (see, e.g., [18] and the references there). Results of Nesterenko [12] show 
that Cg(2) is transcendental for every algebraic number q with < |g| < 1. Zudilin gave 
an upper bound for the measure of irrationality of Cg(2) [22] with 1/g G {2,3,4,...}. 
Furthermore Krattenthaler, Rivoal and Zudilin [1 1] proved that there are infinitely many 
Cg(2n) (and infinitely many Cq{2n +1)) which are irrational, and that at least one of the 
numbers Cg(3), Cg(5), Cg(7), Cg(9), Cg(ll) is irrational whenever 1/g G Z and q ^ ±1. From 
now we only use values of q for which p = \jq G N \ {0, 1} = {2, 3, 4, . . .}. 

In this paper we will show how one can obtain good simultaneous rational approximations 
for Cg(l) ^^^ Cg(2) with a common denominator, which are related to multiple little q- 
Jacobi polynomials and Hermite-Pade approximation techniques. The method that we use 
is an extension of the method that Van Assche [18] used to prove the irrationality of Ci3(l) 
and is an application of Hermite-Pade approximation of a system of Markov functions 
[13, Chapter 4] [19]. 

In Section 2 we will describe the construction of the simultaneous rational approximants 
using Hermite-Pade approximation to two Markov functions which are chosen appropri- 
ately. The solution of this Hermite-Pade approximation problem depends on multiple 
little g-Jacobi polynomials [14], and in Section 3 we give the relevant formulas for the 
special case where all parameters are equal to zero. In Section 4 we will show that these 
rational approximations give a good rational approximation to Cg(l) ^''^^ ^^ particular we 
will prove the following result: 

Theorem 1.1 Suppose p > 1 is an integer and q = 1/p. Let «„ and f3n he given by 
(4-3)-(4-4) for all n G N. Then an, Pn € Z and /9„(^g(2) — «„ 7^ 0. Furthermore 



l/n- 



lim |/3„C,(1) -«„!'" <p-^^ <1. 



In Section 5 we show that we also get a good rational approximation to Cq{2): 

Theorem 1.2 Suppose p > 1 is an integer and q = 1/p. Let an and bn be given by 
(5.2)-(5.3) for all n &'R. Then a„, bn & 'Zj and bnCqi/^) — CLn 7^ 0. Furthermore 



lim \bnCq{2) -an\ ' <P -^ < 1 



These rational approximations are good enough to conclude that Cq{l) and Cq{2) are 
irrational, by using the following elementary lemma: 



Lemma 1.1 Let x be a real number and suppose there exist integers an,bn [n G N) such 
that 

(1) bnX — an j^ for every n eN, 

(2) lim„^oo \bnX - an\ = 0. 

Then x is irrational. 

This lemma expresses the fact that a rational number is approximable to order 1 by 
rational numbers and to no higher order [9, Theorem 186]. The measure of irrationality 
r{x) can be defined as 

r{x) = inf{r G M : \x — a/b\ < l/¥ has at most finitely many integer solutions (a, 6)}. 

It is known that if \x - an/bn\ = O^l/bn^"") with < s < 1 and 6„ < 6„+i < &^+°^^\ 
then the measure of irrationality r{x) satisfies 2 < r{x) < 1 + 1/s (see, e.g., [4, exercise 3 
on p. 376] for the upper bound; the lower bound follows since every irrational number is 
approximable to order 2 [9, Theorem 187]). 

Our rational approximations to Cq{l) give the upper bound 5.0444 . . . for the measure of 
irrationality of Cg(l)) which is not as good as the upper bound of 2.5082 ... in [18] or the 
upper bound of 2.4234. . . in [23]. For Cg(2) we obtain the upper bound 15.8369. . . for 
the measure of irrationality, which is not as good as the upper bound 4.07869374 . . . that 
Zudilin obtained in [22]. Our rational approximations, however, use the same denominator 
and hence we have constructed simultaneous rational approximants. The price we pay for 
this is that the order of approximation for each individual number is not as good as 
possible. But we gain some very important information because our simultaneous rational 
approximants are good enough to prove our main result: 

Theorem 1.3 The numbers 1, Cq{l) and Cq{2) are linearly independent over Q. 

This result is stronger than the statement that Cq{l) and Cq{2) are irrational. We prove 
this theorem by means of the following lemma, which extends Lemma 1.1 (see [10, Lemma 
2.1] for a similar lemma but with a different first condition). 

Lemma 1.2 Let x, y be real numbers. Suppose that for all (a, b, c) E 1j^\ (0, 0, 0) and for 
infinity many positive integers n G A C N, there exist integers p^, Qn (ind r„ such that the 
following three conditions are satisfied: 

(1) apn + bqn + CTn 7^ for every n G A, 

(2) \pnX — g„| — » as n ^ oo, 

(3) IPnV -r„\ ^0 as n^ oo. 

Then l,x,y are linearly independent over Q. 

Proof. Suppose that l,x and y are linearly dependent over Q. Then there exist integers 



oi, 0-2, bi, 62, ci, C2 such that 

ai bi ci 

h —X H y = 0. 

a2 O2 C2 

When we muhiply this by a2b2C2 we get aib2C2 + a2biC2X + a2^2Cit/ = 0, hence there exist 
integers a, b, c such that a + bx+cy = 0. Multiplying by p„ and adding the terms —bqn — crn 
on both sides, we obtain 

On the right hand side we have an integer different from zero, hence the right hand side 
is in absolute value at least 1 for every n G A. But the expression on the left hand side 
tends to as n tends to infinity, hence we have a contradiction. □ 

If we take x = Cq{l) and y = Cq{2), then Condition (2) follows from Theorem 1.1 and 
Condition (3) follows from Theorem 1.2. In Section 6 we will show that Condition (1) of 
Lemma 1.2 holds. This requires some results about cyclotomic polynomials and cyclotomic 
numbers from number theory. 



2 Multiple little g-Jacobi polynomials 

2.1 Multiple little q-Jacobi polynomials 



Little g-Jacobi polynomials are orthogonal polynomials on the exponential lattice {q^, k = 
0, 1, 2, . . .}, where < g < 1. In order to express the orthogonality relations, we will use 
the g-integral 

/ f{x)dqX = J2q'f{q'), (2.1) 

(see, e.g., [1, §10.1], [8, §1.11]) where / is a function on [0,1] which is continuous at 0. 
The orthogonality is given by 



where 



1 

Pn{x;a, l3\q)x'^w{x;a, (3\q) dgX = 0, k = 0,1, . . . ,n — 1, (2.2) 





w{x-a,b\q) = j^^^^^x-. (2.3) 



We have used the notation 

n— 1 00 

(a; q)n = n (1 - «?^)' («5 9)00 = n (1 - «^^)- 

k=0 fe=0 

In order that the g-integral of w is finite, we need to impose the restrictions a,f3 > 
— 1. The orthogonality conditions (2.2) determine the polynomials Pn{x;a,f3\q) up to a 



multiplicative factor. When we define p„(x; a, P\q) as monic polynomials, they are uniquely 
determined by the orthogonality conditions. 

The Rodrigues formula for p„(a;; a, P\q) is given by 



/ „a+n ^a+n—l\n 

(g"+^"+^g)nw(a:,a,/5|g)' 



Pn{x\ a, P\q) = ,„^/^^„^i.„,\..,^ „ ^,^^ D;-iw(a;, a + n, p + n\q) (2.4) 



and an explicit formula by 



„-n „n+a+/3+l 

y 1 H 



,,,(r,aM,)-i-xrf'--^y^j^Q^.^.^ ^^^ 



q; qx 



1 Nn n(n-l)/2 19' i Q)n y^ 



(2.5) 



We have introduced multiple little g-Jacobi polynomials in [14]. Suppose that f3 > —1 and 
that tti, 0:25 • • • , ttr are such that each ai > —1 and ai — aj ^ Z whenever i 7^ j. Then the 
multiple little g-Jacobi polynomial Ph{x; a, P\q) for the multi- index n = (ni, n2, . . . ,nr) G 
N*" is the monic polynomial of degree \n\ = ni+n2 + . . .+nr that satisfies the orthogonality 
conditions 

1 

Ph{x; d, P\q)x''w{x; aj, P\q) dqX = 0, A; = 0, 1, . . . , Uj — 1, j = 1, 2, . . . , r, (2.6) 

These are the multiple little g-Jacobi polynomials of the first kind. We only consider the 
case r = 2, so we take n = {n,m), with m < n. From [14] we know that the Rodrigues 
formula for Pn,m{x; (ai, 0:2), /9|g) is given by 

Pn,m{x; (01,02), f3\q) - 






X D" , [x"^-"2+"D™ iw;(x, a2 + m,l3 + n + m\q)] (2.7) 



and that an explicit expression for Pn,m{x; (cti, 0:2), /9|g) is given by 
Pn,m{x; (ai, ^2), /9|g) = (-1) 



„+^g""^(g"i+^ g)n(g"^+^ g)^g"("-l)/2g"^(>»-l)/2 



.. y- y- (g""; g)j(g"'"; 9)fc(g"^+"+^ g)fc(g"^+^+"+^ g),+fc(g"^+^+"+'"+^ g)^ 
"" to.to (g"^+'; g),+fc(g"^+'3+"+^ g)fc(g"^+^ q)k 

X fej . . ^. (2.8) 

r (g;g)fc(g;g)i 



2.2 The case a\ = a2 



When we take ai = a2 = «, then (2.6) gives only n conditions. However, when we look 
at expression (2.8) we see that Pn,m{x; {ai,a2), P\q) still has degree n + m. We use the 
notation pl^^m'^^x) for Pn,m{x; {ai,a2), P\q) with ai = a2 = a. We will now show that 
PrTm'^'*!^)' defined by the Rodrigues formula (2.7), again satisfies n + m orthogonality 
conditions, namely 



p, 



^^,;:^^Hx)w{x, a, P\q)x'd,x = 0, £ = 0, 1, . . . , n - 1, 



p1^j^''^\x)w{x, a, (3\q)x^ logq X dqX = 0, i = 0,1, . . . ,m — 1. 

To prove these orthogonality conditions we will need the following lemma. 
Lemma 2.1 (Summation by parts) When g{q^^) = we have that 

oo oo 

E //(?') D,-.g{x)\^, = -q^l'diq') DJix)\^, . 

k=0 k=0 



(2.9) 
(2.10) 



(2.11) 



We start with the first n orthogonality conditions. When we use the Rodrigues formula 
(2.7), we obtain 



-1 
/o 
with 



pi';':;^'^\x)w{x, a, l3\q)x'^ dgX = C, 







x"'D'^-iw{x,a + m,[3 + n + m\q) x^ dqX, 



G„ 



f „a+n+m „a+n+'m—l\n f „a+m „a+m—l\m 



(^^/3+a+n+m+l. ^-^^^^p+a+n+m+l.^ ^)^ 



Now we can rewrite the right hand side as an infinite sum using (2.1) and then apply 
summation by parts (2.11) n times. This gives 



Pn,^''^\x)w{x,a,f3\q)xUqX 

oo 

= (-l)"g"a,m E ^' k"D™ iti;(a;, a + m,(3 + n + m\q) 



fc=0 



D'^x^ 



x=q'' 



x=q'' 



Since D^x^ = for i < n, we see that (2.9) holds. 

Now we prove (2.10). Again using the Rodrigues formula (2.7) and applying summation 
by parts (2.11) n times 



Pn,;^''^\x)w{x,a,p\q)xHoggXdgX 

= (-l)"g"C„,m / D"^ (^x^ logg a;) x"'D'^-iw{x, a + m, p + n + m\q) dqX. 



To calculate D" fx^log xj we can use the g-analogue of Leibniz' rule 



D'i[f{;i)g(x)\ 



fe=0 L"'J q 



D;~'f]{q'x)(D'^g]{x), n = 0,1,2,.... (2.12) 



This gives 



piX'^^ (x)w(x, a, P\q)x^ log„ X dgX 



I n 



cl E 



k=0 



(^ 



n-k e 
q -^ 



Dq^Og^X 



x"'D''I'-iw{x,a + m, P + n + m\q) dgX, 



with C a constant. Obviously we have for k < n — £ that V^^ x = and ior k > n — i 
that V^^^x^ is a polynomial of degree £ — n + /c, say nk/-n+k- It is also easy to see that 
Dg log X is a constant times l/a:'^ for /c > 1. So we get 



Pi^m '^Ha;)^^(a^. "> /^k)a;^ logg a: cigX 



„1 n 
/ E Cfc 



nk/-n+k{,q^x)—rx"'D'^-iw{x, a + m,f3 + n + m\q) dgX, 



X 



where the Cfc are constants. Now 'iTk/-.n+k{(l'^x)x"'/x^ is a polynomial of degree i, say Pk/{x). 
Using the Rodrigues formula (2.4) for little g-Jacobi polynomials, we get 



Pim^\x)w{x, a, p\q)x^ log^ x dqX 



k=n~l 



Pfc/(a;)w(x, a, /3 + n|g)p^''^+"^(a;) rf^^x. 



where the c'^ are constants. The orthogonality relations of the little g-Jacobi polynomial 
imply that the right hand side is for i < m, so that (2.10) holds. 



2.3 The case cti = 0:2 = /3 = 



From now, we will work with the multiple little g-Jacobi polynomial Pn,m{x', («!, 0:2), /9|g) 
where ai = a2 = f3 = and with another normalisation, namely 



Pn,m(0; {ai,a2),P\q) = 1. 



(2.13) 



We will denote this polynomial by Pn,m{x). Since tf;(a;, 0, 0|g) = 1, it follows from (2.9)- 
(2.10) that pn,m can be defined by 



Pn,m{x)x^dqX = 0, i = 0,1, ... ,11 - 1, 

Pn,m{x)x^ logg xdqX = 0, £ = 0, 1, . . . , 771 - 1, 

together with the normahsation (2.13). 
From (2.8) we see that 



(2.14) 
(2.15) 



Pn,m(0; {ai,a2),(3\q) = (-1) 



n+mg""!?"^^'; g)n(g"^+^ q)^qn(n-l)/2qm(m~l)/2 



(q^ 



+ai+n+m+l ■ 



qUq^ 



+a2+n+m+l 



;q)r 



(2.16) 



Multiplying (2.7) and (2.8) by [p„,m(0; {ai,a2), P\q)] ^ and setting ai = a2 = P = gives 



and the Rodrigues formula becomes 



PuA^) = (-l)"+-g"("-i)/^g-(— 1)/^ ^^ /_ D;-.[x"D-x ((gx;g)„+„a:^ 

W) Q)n{q-j q)m 



f2.18) 



Equation (2.17) expresses the polynomial p„_m in the basis {1, x,x^, . . . , x"-^"^}. Sometimes 
it is more convenient to use the basis {{qx; q)^, < i < n + m}. The Rodrigues formula 
(2.18) allows us to obtain an expression for Pn,m{x) in this basis. Recall the g-analog of 
Newton's binomial formula 



{x;q)n = Yl 

fc=0 



gk(k^l)/2^_^-^k^ 



(2.19) 



and its dual 



fc=0 



(_l)fcg-«fc+fe(fe+l)/2(3,.^)^_ 



This is a special case of the g-binomial series [1, §10.2] [8, §1.3] 

|g| < 1, \x\ < 1. 



y, (a;g)n ^n _ {ax;q)c 



Using (2.20) with argument q'^~^"^~^^x and exponent m gives 



(2.20) 



(2.21) 



j.m _ (n+m+l)m y^ 
fc=0 



(-l)kq~mk+k{k+l)/2^^n+m+l^. ^^^ 



Using this in the Rodrigues formula (2.18), we find 



(n — 1 )n+m 
( \ / -1 \n+m n(n—l) /2 m(m—l) /2 V^ / (n+m+l)m 



iQ;q)niq;q), 



.py 



E 

fc=0 



m 
k 



'^^l)kq-mk+kik+l)/2j^n_^^^njjrn_^(^^^.^^ 



n+m+fcj ) 



because {qx; q)n+m (g"'"'"™^^x; q)k = {qx; q)n+m+k- One easily finds that 



Dp{qx;q), 



{q;q)n-k{l-p)^ 



{qx; q) 



n—ki 



f2.22) 



so we get 



(n — ^ "1"+™ 

„ (^\ — ( 1 \n+m n{n-l)/2 m(m-l)/2 K'd ^J (n+m+l)m 



(y;y)n(g;g), 



.pv 



xE 

fc=0 



m 
A; 



-I? 



/ iNfc -mfc+fc(fc+l)/2 (g;g)n+m+fc n« \^"■(nT■n^ ,1 



Using (2.20) again, this time with argument q^^^ x and exponent n, gives 



x" = p("+^+i)" y^ 



j=0 L^Jg 



.l)i^-nj+ja+l)/2(^n+fc+l^.^)^.^ 



and since (gx; q)n+k (g"^ x; g)j = (gx; q)n+k+j we obtain 



Pn,m\-^ ) 



1 \n+m n(n-l)/2 m(m-l)/: 



l)/2 W -"-J „(n+m+l)m V^ 



(y;y)n(g;g), 



.pv 



772 



fc=0 L^Jq 



-1)^9 



fc -mA:+A:(A;+l)/2 



(y;y)r 



(y;g)n+fc(i-p)™ ,t^o 



(_l).g-n.-H-.-a+l)/2nn_^(^^.^)^^^^^.. 



Again using formula (2.22) and also using that 



{q;q)k = i-l)'p~'^'^'^^'{p;p)k, 



n 
k 

L Jq 



P 



-k(n—k) 



we can rewrite the expression for Pn,m{x) as 



Pn,^(x) = (-1)"+™EE(-1)'^^' 
k=Oj=0 



n + m + k 
m 



n + k + j 



n 



- p 


m 
k 


p 


n 



X (ga;;g)fc+,p(""^')("-^'+^)/V"^^'=)('"-'=+^)/^ (2.23) 



3 Hermite-Pade approximation of /i and /2 



We define two measures /ii and fi2 by taking dfii{x) = dqX and dfi2{x) = \ogg{x)dqX, 
wliere dg is defined by (2.1). Tlien /ii and fi2 are supported on {q'', k = 0,1,2, . . .}, wliicli 
is a bounded set in [0, 1] witli one accumulation point at 0. Tlie Markov functions for tlie 
measures /ii and /i2 are 



/•I dfii{x) 

/ Z — X 


^ q' 
^ z - a'' 


/•I dfi2{x) 

z — X 





M^) = L -^^ = E — ^- (3.2) 

Observe tliat for every N &N 



ui)=W)+i: 



N-l 1 



fc=lP''-l' 



^-1 A; 



C,(2) = Mp^)+E-nrT + ^^i^^' 



therefore we will look for rational approximants of fi{z) and f2{z) with common denom- 
inator and evaluate these at z = p^ for an appropriate choice of N . We can find these 
rational approximants by Hermite-Pade approximation of type II. 

For Hermite-Pade approximation of type II one requires a polynomial Pn,m of degree 
<n + m, and polynomials g„^m and r„^m such that 

Pn,m{z)fi{z) - qn,m{z) = O f ^^^j , Z^OO, (3.3) 

Pn,m{z)f2{z) - rn,m{z) = O f ^^;q:Y j • Z ^OO, (3.4) 

It is known [13, Chapter 4] that for m < n the polynomial Pn,m, is, up to a multiplicative 
factor, uniquely given by 

1 

Pn,rn{x)x'' d^li{x)=0, £ = 0, 1, . . . , U - 1, (3.5) 



1 

Pn,rn{x)x'^ d^2{x)=0, £ = 0, 1, . . . , m - 1, (3.6) 

/O 

and that qn^m and r„ ,„ are given by 

qn,m{z) = dfli{x), rn,m[z) = I dH2{x). 

Jo z — X Jo z — X 

(3.7) 



10 



The remainder in the approximation (3.3) is 



Pn,m{z)fl{z) -qn,m{z) = I -^^^ d^i{x), 



Z — X 



(3i 



and for (3.4) 



Pn,m{z)f2{z) -rn,m{z) = I -^^^ dll2{,x). 



Z — X 



(3.9) 



Comparing (3.5)-(3.6) with (2.14)-(2.15) we see that the common denominator is given 
by the multiple little g-Jacobi polynomial Pn,m- So (2.23) gives an explicit expression for 
Pn,m- Then it follows that we can compute qn^m and r„,m explicitly using (3.7). For g„m 
we have to compute 



qn,m{z) = 2^ 2 q 

By using the explicit expression (2.23) for Pn,m{z), we find 



(3.10) 



m n 



g.,^(z) = (-l)"+™EE(-l 

k=Oj=0 



n + m + k 
m 



n + k + j 



m 



k 

Jp L jp 



Jp 



X p 



n 

-ip 

{n-j){n-j+l)/2{m~k){m-k+l)/2 Y^ (O'^J Q.)k+j — [Q ] Qjk+j ^l 



p. 






z — q^ 



Now use 



^''''^'^'-^'''^'^ = -t,\w,,).M-^x-.,)k^. 



which one can prove by induction, then this gives 



qn,rn\Zj 



m n 
\n+m+l \"^ \'^(' 1 \k+j 



EE(-i: 

A;=Oj=0 



n + m + k 
m 



n + k + j 



n 



- p 


m 
_k_ 


p 


n 
J. 



k+j 



X p{n~j){n~j+i)/2p(m-k)im-k+i)/2 J- q^ {q^+h; q),„,^^ ^ q'{q'+'; q),^,. 



r=l 



By using the g-binomial series (2.21), we can compute the modified moments 






(3.ii; 



so that 



m n 
n+m+1 \ ' \ '/' -| \k+j 



gn,™(^) = (-i)"+™+^EE(-i 

X p 



n + m + k 
m 



n + k + j 



n 



m 
k 



k+j / f+1 . A 
(n-j){n-j+l)/2(m-k){m-k+l)/2 V^ W Z] q)k-r+j /g ^^^ 



11 



For an explicit expression of rn,m we use (3.7), which gives 

Pn,m\Z) Pn,rn\C[ ) /, £ 



I n.in, \ 2- ) 



E 



z — q" 



(3.13) 



Completely analogous to g„^m we find 



' n.m\ ■^ ) 



m n 
\n+m+l \ ^ \ ^( -| \k+j 



EE(-i; 

fc=oi=o 



n + m + k 
m 



n + k + j 



k+j 



n 



jp 



m 
k 



n 



? 



r=l ^=0 



Now we have more work to compute 



oo oo / r. „\ 

EV(?^-^^?).-i = (?;?).-! EVf^ 



Using the g-binomial series (2.21), we have that 



£=0 Wi^J^ 



By taking the derivative with respect to x we find 



d_ 
dx 



{<f\ q) 



[(a;; <i)r] E / \ ^^ + (2^; ?)r E 



£=0 



(g;g) 






„^-l 



x^-^ = 0. 



It is not difficult to see that 



so we find 



5^1'"' "'■ 



ff(^x' = . 



Eyx, q)r 
1 — i " 



j=0 



xg' 



r-l 



E— ^E 



~ (g-g) 



% ii; Q)i ^ 1 - a;g^ f^Q (g; g)^ 

From this it follows that (3.14) becomes 



X . 



r-\ 



Y^iq^q^^^; q) r-l = {q;q)r-iqJ2 T 

n .■ — n -L 



g 






g^+^ ^ (g; q)i 



£=0 i=0 -" 'i £=0 

Using expression (3.11) for the modified moments we find 



QO 1 *" 1 



g' n P' 



(3.14) 
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When we use this in the expression for rn,m we get 



fc=Ojr=0 



n + m + k 
m 



JpL 



n + k + i 



n 



m 



JpL JpL^Jp 



n 



k+j r ( r+1 . \ 



fl^(P'--l)(P*-l)' 



(3.15) 



We will evaluate these functions Pn,miz), qn,m{z), rn,m{,z) dX z = p"-+'". First we will show 
that Pn,m{p"'~^"^) is an integer. From the g-version of Pascal's triangle identity 



n 




n — 


ll 


+ / 


[n-ll 




[n-ll 


k 


- 


k- 


1 


k 


- 


k 




p 






p 




p 





+ p 



n—k 



n — I 
k-1 



(3.16) 



jp 



it follows (by induction) that ? is an integer whenever p is an integer. Furthermore we 



have 



(gp"+"^;g)fe+, = (p 






i=l 



For each value of i in this product, the factor 1 — p"+™ * is an integer. This means that 
Pn,m(p"+"'), with pn,m giveu by (2.23), is an integer. 

Now we will evaluate qn,m{z) and r„_m(z) at p"+'". We can use (3.12) at z = p"-+™ and 

(g"+ V+"^; g)fc-r+, = (p"+™-^-^;p)fc-^+„ to find 



Hn,m\P 



ra+m\ 



m n 



EE(-i: 

A;=Oj=0 



n + m + k 
m 



-ip 



n + k + j 
n 



m 
k 

Jp L jp 



Jp 



^+3 (^n+m-k-j. \ 
X p(n-i)(n-i+l)/2 {m-fc){m-fe+l)/2 V- U^ ,P)k-r+j ^ (3.17) 

r=l p*" — 1 

The terms in the sum for qn,m{p"'^"^) are not all integers, because of the expression p^ — 1 
in the denominators. In order to obtain an integer we have to multiply qn,m{p"'~^"^) by a 
multiple of all p"^ — 1 for r = 1, 2, . . . , n + m. 

For rn,mip"'~^™') we use (3.15) at z = p"+™' and again (^^'■+1^"+™- g)^_^_^^. = (p"+"*~'^"-5''; p)^.^.,.^-, 
to find 



'n,m\P 



n+rn\ 



k=Oj=0 



n + m + k 
m 



n + k + j 
n 

k+j r ( n+m-k-j. 



- p 


m 
k_ 


p 


n 



X p{n~j){n-j+l)/2 (m-k){m^k+l)/2 V- y- KP j_P)fc-r+i ^ (3.18) 



{p 



We see that in order to get an integer now, we have to multiply rn,m{p^~^^) by a multiple 
of {p"^ — l)(p* — 1) for r = 1,2, . . . ,n + m and i = 1,2, . . . ,n + m. 
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Define 

d^{x) = f[^kix), (3.19) 

fc=i 

where 

n 

$„(a;) = n (^ - e'"""^") (3.20) 

fe=i 

gcd(fc,n)=l 

are the cyclotomic polynomials. Each cyclotomic polynomial is monic and has integer 
coefficients. It is known [17] that 

x''-l = llM^), (3.21) 

d\n 

and that every cyclotomic polynomial is irreducible over Q[x]. Hence dn{p) is a multiple of 
all p^ — 1 for i = 1,2, ... ,n. The growth of this sequence is given by the following lemma 

[18, Lemma 2]. 

Lemma 3.1 Suppose x > 1 and let dn{x) he given by (3.19). Then 

lim dJxY/''' =x^'^\ 



Since dn+m{p) is a multiple of p^ — 1 for all £ = 1,2, ...,n + m, we conclude that 
c^n+m(p)gn,m(p"^™) and dl^^{p)rn,m{p"'^"') are integers. 



4 Rational approximations to Cq{l) 



From (1.1) we know that 



oo n 






In this section we will construct rational approximations to Cq{l), prove the irrationality 
and give an upper bound for its measure of irrationality. 



4.1 Rational approximations 

In (3.1) we defined the function /i. When we evaluate /i at p"+™', we get 

OD k 00-1 n+m—1 1 
J l\y ) Z_^ ^n+m ^k ^-^ ^n+m+k 1 ^1^ ' 2-^ „fc 1 ^ ' ' 

k=oy y fc=o^ -^ k=i y ^ 
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and hence 



n+m— 1 



C,(i) = /i(p"+'") + E 



k=l 



pfe _ X 



We now take m = n — 1 and define 



(3n = d2n-l{p)Pn,n-l{p " ) 



an = d2n-l{p) 



2n-2 



gn,n-l(/"-')+Pn,n-l(/"-') E 



fc=l P'' - 1, 



(4.2) 
(4.3) 

(4.4) 



Then it follows from equation (4.2) and from the Hermite-Pade approximation of /i (3.8) 
that 



/5„Cg(l) - an = d2n-lip)[Pn,n-liP " )/l(P " ) - qn,n~lip 

, / X /"^ Pn,n-l{x) 

= d2n-l{p) j^ ^^=r^dgX. 



2n-l\ 



(4.5) 



4-2 Irrationality of (q{l) 



Using Lemma 1.1, we can prove the irrationality oi (q{l). We will first show that 

lim IPnCgC^) -«n| = 0, 

i.e., using (4.5), we want to find an estimate for 

[^ Pn,n-l{x) 

/ — ■ d„x. 

Jo z — X 

When we use the Rodrigues formula (2.18) for pn,m we get 



where 



^^^^^ dqx = C f' D^-,[x-D^., {{qx- g)„+„x'")]^ dqx, 
z — X Jo ^ ^ Z — X 



(n — 1 ')"+'" 
/^ / -i\n+m n{n—l)/2m{m—l)/2 Vi -^ J 



{q; <l)n{q; q) 
Repeated application of summation by parts (2.11) gives 



Pn,m\X J 



10 Z ~ X JO 

Now it is easy to see by induction that 



Dl 



z — X 



dgx = C{-qr / a:"D"l 1 {{qx- g)„+™x™) Dl 



{q;q)r 



z — X 






(g;g)„(l-g)' 



{z — x){z — xq) . . .{z — xq^) (1 — g)"2;"+i (-; q) 



n+l 



(4.6) 
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so for z = p"+'" we find 



p"'+™ — X 



dgx = C(-g)" 



(g;g)r 



(1 - q)np{n+m)in+l) J^ {q^+^X] q) 



x"D™i {{qx\q)n+mX'" 



n+l 






Again we apply summation by parts (2.11) several times and use the g-analogue of Leibniz' 
formula (2.12). This gives 



Pn,m\X) 



p 



n+m 



X 



dqX = C{—q) 



n+m 



{q\<i)r 



'\ Q\nji(n+Tn){n+1) 



m 



(ga;;g)„+„x™D™-'=(x") 



k=0 L ^q' 

By induction, it is easy to see that for r < s 



D 



k^ « (g^+^x; g)„+i ^' 






Dl {x') 



iq;q)sX' 



and that 



Using this we find 



D 



[l-qYiq;q)s-r' 

g'=("+'")(g"+i;g)fe 



(4.7) 



•? (g^+'^x; g)„+i (1 - g)'=(g«+-x; q)n+k+i 



Pn,m[X) ^ _ (^(_„\n+m I?' 9)n 



m 

xE 



m 



k=0 L^Jg 



Q 



k{2n+k) 



(9;g)n(g""^^g)fe r^ {qX;q)n+mX 



n+k 



dqX. (4.8) 



'1 - ?)""(?; q)n-m+k Jo (g^+^x; g)„+fc+i 
The integrand is always positive and (^"'^'"x; q)n+k+i > (g^+^j g)„_,_fc_,_i, hence we find 



(gx; g)„+mX"+'= 



(i„x < 



(g"+™x; q)n+k+i ' ~ (g"+'"; g)n+fc+i Jo 



yqXj qjn-^mX (XqX 



< 



< 



(g"+™; q)n+k+i JZq 



E(^^^^'^)«+'»(^ 



t\n+k+l 



[Q'lQjn+m Y> (g" ™ 'iQJe ,n+k+i\e 



(g"+™; g)„+fc+i ^ (g; q)e 
Using the g-binomial series (2.21), we then find 



1 (gx; g)„+^x"+'= ^ ^ (g; g) 



/'„2n+fc+m+2. \ 
n+m Vy ) "i/c 



(g"+-x; g)„+fc+i " (g«+-; g)„+fc+i (g"+'=+i;g)c 



(g;?) 



n+k 



(rin+m- ri^ , _ (rin+m+1- ri^ , _ ' 

[Q jQJn+k+iyq ^ ^ ,q)n+k+i 
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This estimate, together with the definition of C in (4.6), then gives 



Pn,m\-'^ ) 



p 



:n+m 



lA/fltAj 



X 



< 



„n{n—l) /2 „m{ra—l) /2 „n+m 



m 



fe=0 l^q 






+fc+i 



Some simple estimations then give 



Pn,m\X ) 



P' 



Xi+m 






X 



< 



n{n~l) /2 ra(m—l) /2 n+m 



,2?i 



pin+m)in+l) (^. ^)^ (^n+m. g)2_^^^2 fro ^ 
„ra(n— l)/2„m(m— l)/2„n+m i „2?i(?Tt+l) 



2n 



This gives a useful estimate for the integral on the right hand side of equation (4.5), which 
(for m = n — 1) implies that 



\PnCq{l) -«n| < d2n-l{p) 



^n(n-l)/2^(n-l){n-2)/2 



,2n-l 



p(2n-l)(n+l) (g. g)„_,(g2n-l. ^)2^^^(1 _ ^2n) ' 



(4.9) 



Now we can prove Theorem 1.1 

Proof of Theorem 1.1. From Section 3 we know that «„ and /5„ are integers. From 
the equations (4.5) and (4.8) it follows that /9„Cg(l) ~ ctn 7^ because d2n-i{p) 7^ and 
because the integral on the right hand side of (4.5) can be written as a sum with all terms 
different from zero and of the same sign. From (4.9) we can find that 



|l/n 



lim \(3r.Ul) - a^r' < lim d^n-iipr'' q 



l/n^ 3 



..p^V-'p-^ 



p 



<1, 



where we have used Lemma 3.1. Hence lim„^oo IPnCqi^) 
follows from Lemma 1.1. 



— ar 



0. The irrationality now 



D 



4-3 Measure of irrationality for (g(l) 



Theorem 1.1 gives rational approximations an/Pn for C<?(1) that satisfy 



ai)-f^ 

Pn 



o 



^^(-^^"^\ 



Pn 



(4.10) 
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for every e > 0, with /?„ = d2n-i{p)Pn,n-i{p'^"'~^) ■ We already know the asymptotic be- 
havior of d2n~i{p), so what remains is to find the asymptotic behavior of Pn,n~i{p'^"'~^) as 

n — s> oo. 

Since {l,x,x^, . . . ,x"'~'^,\ogx,x\ogx,x'^\ogx, . . . ,x"^~^\ogx} is a Chebyshev system on 
[0, 1] whenever m < n, we know that all the zeros of Pn,m{z) are simple and lie in the 
interval (0, 1). If we call these zeros xi, X2, . . . , Xn+m then we can write Pn,m.{,x) as 

n+m 
Pn,rn\X ) rin,m J_J_ \-^ Xjj. 

i=i 
For |a;| > 1 we have |x| — 1 < |a; — Xj\ < \x\ + 1, hence 

n+m 



{\x\ - 1)"+™ < J] \X-Xj\< {\x\ + 1) 



n+m 



i=i 



so multiplying by \nn,m\ and evaluating at x = p""*"™ gives 
For m = n — 1 one obtains 

2ra— 1 2n-] 

SO we have that 



lim \pn,n-i{p 



,2n-l\|l/n2 



p*^ lim I /€„,„-! 



l/n2 



Equation (2.16) implies that the leading coefficient of j9„^m is given by 



K, 



n,m LFn,m 



.o,o)/n^l-i 



-1) 



n+m 



q''"'{q; q)n{q] g)„g'^('^-i)/2gM™-i)/2 • 



(4.11) 



Since {q^ , q)k = (q', q)n+k/iq', q)n, we can rewrite the leading coefficient of p„,m as 



Kr, 



\n+m 



This gives lim^^oo |/t„,„_i 



l/n^ 



2n + m 




n + 2m 


n 


9 


m 


p^, so we 


h 


ave 



(n+m)V2 -(n+m)/2^ 



lim K„_i(p'"-^)r/"'=/. 



Combining this with Lemma 3.1 we have for the integers /3„ in (4.3) that 



6(7r"^+2) 



1- I/O |l/n" 12/7r2 6 2iZL+ii 

lim |p„| ' = p ' p = p^ . 



(4.12) 



(4.13) 



Together with (4.10) this gives that 



C,(i) - f 

Pn 



O I -—^^^^ I (4-14) 



for every e > 0, which imphes the following upper bound for the measure of irrationality 

.(C,(l)) < 1 . '-K^ - 4^ « 5.04443. 

TT^ — 4 TT^ — 4 



5 Rational approximations to Cg(2) 
From (1.1) we know that 

oo ^^n oo ^ 



ra=l 



g" ^=1 P" 



In this section we will construct rational approximations to Cg(2), we will prove that Cq{'^) 
is irrational, and we will give an upper bound for its measure of irrationality. 



5.1 Rational approximations 



Evaluating the function /2, which we defined in (3.2), at the point p"+™, gives 



oo I. k oo 7, 

f J n+m\ ^ y- '^H ^ y- "^ 



k=0 f ^ fc=0 

So we have 

/2(p"+-) = Vi^±^^^-(n + m)E— -4i 

A;=0 " fc=0 " 

n+m— 1 7, 

= C.(2)- E ^-r-{n + m)h{p-^n. 
where we used (4.1) for /i(j5"+™). Hence we can write Cg(2) as follows: 

n+m— I 7, 

C,(2) = /2(p"+'") + E ^3- + (n + m)/i(p"+-). (5.i; 

k=i P -"- 

We now take m = n — 1 and define 
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&n = rfL-l(pKn-lb'"-'), (5.2) 

2n-2 7, 



a„ — «2n-l(P) 



fc=l P - 1 



. (5.3) 



Then it follows from (5.1) and from the Hermite-Pade approximation of /i and /2 (3.8)- 
(3.9) that 



6„Cg(2) -an = rf2„_i(p) 



-^^^^^ rf/i2(x) + (2n - 1) / ^^3^ ^ rf/ii(a;) 



P 



X 



p. 



X 



(5.4) 



From Section 3 we know that the numbers a„ and 6„ are integers for all n G N. We will 
show that 

lim |6„Cg(2) -a„| = 0, 

and that 6„Cq(2) — a„ 7^ for all n G N, so that Lemma 1.1 implies the irrationality of 
Cg(2). But before we can do that, we have to prove some results about the asymptotic 
behavior of pn,m as n, m ^ 00. 



5.2 Asymptotic behavior of Pn,, 



It is well known that the common denominator of Hermite-Pade approximants satisfies 
a multiple orthogonality relation. Driver and Stahl [5] showed that for a Nikishin system 
these common denominators also satisfy an ordinary orthogonality relation. Although in 
our case (/i, /2) do not form a Nikishin system, we can prove in a similar way that Pn,m in 
our case also satisfies an ordinary orthogonality relation. We need the following theorem. 



Theorem 5.1 Let pn,m be the multiple little q-Jacobi polynomial given by (2.14)-(2.15), 
^^t Qn^m be defined by (3.10) and let m < n. Then we have 



{y - '^f\pn,m{y)fi{y) - qn,ra{y)] dy = o 



(5.5) 



for k = 0,1, . . . ,m — 1. 

Proof. From the expression of the remainder of the Hermite-Pade approximation (3.8) 
it follows that 



(y - ^)''[Pn,miy)fiiy) - qn,miy)] dy 



{y-i 



\k+l 



y-1 



Pn,m\-^j 

y — X 



dfii{x) 



dy. 



If we add and subtract {x — 1)^^^ on the right hand side, then we have 
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(y - ^f[Pn,miy)fiiy) - qn,miy)] dy 







+ 



ooy-1 
1 



{y - l)'=+i -{x- 1)^+1 

y — X 

1 {x - l)'=+i 



■Pn,m{x) dlXi{x) 



dy 



CO 



y-1 



y-x 



Pn,m{x) dl^l{x) 



dy. 



Now Uy— 1)^^^ — (a; — l)^^^j /(y — x) is a polynomial of degree A; in x, so because of 
the orthogonality (2.14) the first integral on the right hand side vanishes for k < n. Since 
m < n, this integral therefore vanishes for k = 0,1, . . . ,m — 1. Changing the order of 
integration gives 



rl rO (x — 1) 

{y-lf[Pn,m{y)fl{y)-(ln,m{y)]dy= {x-lfpn,m{x) 7 --7 T 

Jo J -^ [y - l)[y - x) 

By a partial fraction decomposition we obtain 

1 1 {x-l) 



dyd^i{x). 



1-y x-y {y-l){y-x) 



so we find that 



(x-l) 



dy 



lim 



^ {y-'^){y-x) 

Using this in the previous expression, we find 



dy 
t l-y 



dy 

t X — y 



logx. 



{y - ^)''[Pn,m{y)fl{y) - qn,m{y)] dy= / {X- l)''p„,m{x) logX dfli{x). 



The orthogonality relations (2.15) then imply that the right hand side is equal to zero for 
k < m. □ 

From this theorem it follows that Pn,m{y) fi{y) — (ln,m{y) has at least m sign changes on 
the interval (—00, 0). For suppose Pn,m{y) fi{y) — (ln,m{y) has only r < m sign changes on 

(— cxo, 0), say at si, S2, ■ ■ ■ , Sr, then we have with Tiriz) = {z — si){z — S2) ■ ■ ■ {z — s^) that 



TTrd/) [Pn,miy)fliy) - qn,miy)] dy ^ 



because the integrand has no sign changes on (— oo,0). This gives a contradiction with 
the fact that this integral can be written as a linear combination of integrals of the form 



{y - i)^bn,m(i/)/i(y) - qn,m{y)] dy 



with k < r < m, because from Theorem 5.1 we know that these integrals are all zero. So 
Pn,m{y) fi{y) — qn,m{y) has at least m sign changes on (—00, 0). The condition m <n and 
(3.3) guarantees that these integrals are finite. 
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Let Si, S2, . . . , Sm. be m points where pn,m{z) fi{z) —qn,m{z) changes sign on (—00, 0). Then 
we define the polynomial w by 

w{z) = {z- sx){z - S2) ... (2; - s„). (5.6) 

Now we can prove that j)n,m also satisfies an ordinary orhogonality relation. 

Theorem 5.2 Letw he defined by (5.6) and letpn,m the multiple little q-Jacobi polynomial 
defined by (2.U)-(2.15). Then 



/ Pn,m{y)y 
Jo 



k df^ijy) 

w{y) 



0, 



for k = 0,1, . . . ,n + m — 1. 



Proof. Let 7 be a closed positively oriented path of integration with winding number 1 
for all its interior points, such that the interval [0, 1] is in the interior of 7 and the zeros 
of w are outside 7. Using (3.8) we have that 



1 



y I Pn,m\-^ ) 



djii{x) dy 



1 



27ri Jj w{y) Jo x — y r ^\ j a 2711 j-y 
Changing the order of integration on the left hand side gives 
. . 1 



y^[Pn,rn{y)fi{y) -qn,m{y)] 



Pn,m\^ J 



V dv 1 r 

o • / --d^i{x) = TT- y^[Pn,m{y)fi{y) - qn,r 



dy 

w{y)' 



dy 



/o 'zm Jy X — yw[y) zm j^ ' ' w{y) 

The only singularity inside 7 is a;, hence Cauchy's formula on the left hand side gives 



T 1 r (ill 

Pn,m{x)^r^dfXi{x) = TT- / y''[Pn,m{y)fl{y) - qn,miy)] ^T^- 

w{x) zm J-y u][y) 

Since all zeros of w are also zeros oi pn,m{y) fi{y) - qn,m{y), the function y^[pn,m{,y) h{y) ~ 
qn,m{y)]lw{ii) is analytic in C \ [0, 1]. Furthermore 7 encloses the interval [0, 1], hence for 
R sufficient large 



1 , ^ k dfii{x) 

Pn,rn\X)X ~, r~ 

w[x] 



y^[Pn,miy)fiiy) - qn,miy)] 



dy 



2m Jvn w{y) 

with Tfi the circle with center and radius R. We can estimate the expression by 

kPn,m{y)fliy)-(ln,m{y) 



1 fc dnijx) 

Pn,m\X)X 
W[X] 



< R max 



y 



w{y) 



Using (3.3), we know that the function y^ [Pn,m{y)fi{y) — qn,m{y)]/u!{y) has a zero of order 
n + TTi + l — kat infinity. So when R tends to infinity we find 



1 fc dij,i{x) 

Pn,m\X)X , , 

W\X) 



O 



/?"+ 



m—k 
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Therefore 



1 j^ rf/ii(x) 



10 W{X} 

when k <n + m — 1. □ 

Now we have proved that Pn,m satisfies an ordinary orthogonahty relation. Using this, we 
can easily show that Pn,m{,z) fi{z) — qn,m{z) has exactly m sign changes in the interval 
(— oo, 0), namely at the zeros of w. Using (3.8) we have 

Pn,m[z)jl[z)-qn,m[z) = -—-dfli{y). 

Jo z-y w[y) 
If we add and subtract w{z) on the right hand side, then we get 

Pn,m{z)fl{z) - qn,m{z) = / r~\ dfliiy) + w{z) / -— . 

JO z — y u][y) JO z — y w[y) 

The first integral on the right hand side vanishes because of the orthogonality (Theorem 
5.2), so that 

/ x/./ X / X / .Pn,m{z) [^ Pn,m{y) d^i{y) 

Pn,m[z)ji[z) - qn,m{z) = w{z)- 



Pn,m{z) Jo z-y w{y) 
Now we add and subtract Pn,miy) on the right hand side, so we get 

Pn,m{z)fl{z) ~qn,miz) = 

W{z) f^Pn,miy). ,x_ dfjM_ ^iz) [^ Pi,mi.y) djliiy) 

'\Pn,m\Z) Pn,m\y ) ) / X """ 



Pn,m{z) Jo Z-y w[y) Pn,m[z) Jo Z-y w[y) 

Again, the first expression on the right hand side vanishes because of the orthogonality 
(Theorem 5.2), so we find 

Pn,m{z)fi{z) - qn,m{z) = / -^^^ ]-—. (5.7) 

Pn,m{z) Jo Z-y w{y) 

The integrand is not identically zero and has a constant sign for z e M \ [0, 1], hence 
the integral can not be zero, so it follows from this expression that the sign changes of 
Pn,m{z) fi{z) — qn,m{z) in M \ [0, 1] are at the m simple zeros of w in (—00, 0). 



5.3 Irrationality of (q{2) 



We will now show that lim„_+oo(^nCi3(2) — o-n) = and 6„Cg(2) — a„ 7^ for all n G N, so 
that Lemma 1.1 implies the irrationality of Cq{2). First we can reduce the right hand side 
of (5.4) as follows. When we look at the Cauchy transform oi Pn,m{z) fi{z) — qn,m{z) and 
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use (3.8), then we find 



dz 



Because 



y- z 
1 



(^Z I Pn,m\-'^ ) 



oo y — Z Jo Z — X 



■ djj,i{x) 



{y-z){z-x) (y-x) 
we find by changing the order of integration 



y — z X — z 



Pn,m\Z ) Jl\Z ) Qn,m\Z ) i I Pn,m\X ) , , 

az= I — ■ djii[x) 



y- z 



y — X 



y — z X — z 



dz 



I _n^ \ogx dfii{x) — \ogy / -^^^ dixi{x). (5.8) 

Jo y — X Jo y — X 



By multiplying both sides by 'i^+m(p)/ log Q ^^^ using log x dfii (x)/ log q = log^ x djii {x) 
dfi2{x), we obtain 

dl+m{p) [^ Pn,m{z)fi{z)-qn,n,{z) 



logg 



dz 



y- z 



dn+m\P) 



^Pn,m{x) logy [^Pn,m{x) 

y — X log q Jo y — X 



Evaluating dX z = p"+'^ gives the right hand side of (5.4) when m = n — 1, so it turns out 
that 

6,.C,(2) - .„ = %M I" P....-i(^)/.M-.>...-(^) ,^_, (5,9) 



logg 



p2n-l _ ^ 



Now we will estimate the expression on the right hand side. Multiplying and dividing by 
w{p'^^™') and adding and subtracting w{z) gives that the right hand side is equal to 



dln-i{p) 



logg w{p'^"- ^) 



° / / 2n-U / ..Pn,n-l{z)fi{z) - qn,n-l{z) 

{w{p )-w{z)) ^^ dz 

-oo p — Z 

, /■° ^ •,Pn,n-l{z)fl{z)-qn,n-l{z) , 

+ / w{z) ^^ dz 



The first term is because of the orthogonality relations given by Theorem 5.1, so we 
find 



&nCg(2) - a^ 



dL~iip) 



logg tf;(p2n~i) 



° f ^ Pn,n-l{z)fliz) - qn,n-l{z) ^ , . 

^(^) „2n-l ^ ^^- (5-10) 

-oo p — Z 



Since w is a monic polynomial of degree m = n — 1 and all the zeros of w are in the interval 
(—00,0), we have that w(p^"~^) > (j92n-i^n-i_ ^^ g^jg^ y^q^q i^j^^^t p^n-i — ^ > j9^"~i for 
z G (— oo, 0), so we find 



log g ^ 



Wiz) [pn,n-l{z) fi{z) - qn,n-l{z)] dz 



(5.11) 
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The integral on the right hand side can be evaluated exactly, which is done in the following 
lemma. 

Lemma 5.1 Suppose that m < n — 1, then 

f w{z)[p^,U^)Mz) - qnA^)]dz = g"(-+^)(g; q^ ^HitW ^ I'" l^S^- (5-12) 

Proof. Multiply both sides of equation (5.8) by w{y), then add and subtract w{z) on 
both sides of this expression and use the orthogonality relations (5.5) for the left hand 
side and (2.14)-(2.15) for the right hand side, to find 



[Pn,m{z)fi{z) - qn,m{z)] dz 

-oo y — z 

1 w{z) 



Pn,m{z)\ogzdfli{z) -logy / Pn,m{z)dfli{z). 

Jo y — z 



y-z Jo y 

Multiply both sides by y and add and subtract z in the second term on the right hand 
side to find 

w{z)[pn,„Xz)fl{z) - qn,miz)] dz = / w{z)pn,m{z) \og Z dfli{z) 

-oo y — z Jo y — z 

-logy / W{z)pn,m{z) d^i{z) - logy / w{z)Pn,m{.z) djli^z). 

Jo Jo y — z 

We know that tu is a polynomial of degree m. For m < n — 1 the orthogonality relations 
(2.14)-(2.15) imply that the second integral on the right hand side is 0. Taking the limit 
for y going to oo, we get 

rl 

Wiz)[pn,raiz)fi{z) - qn,raiz)] dz = w{z)p.a^rniz) \og Z dlii{z). 

-oo Jo 

Since w is a monic polynomial of degree m, it follows from the orthogonality relations 

(2.15) that 

rl 

W{z)[pn,m{z)fliz) ~qn,ra{z)]dz= z"" p n, m{z) \og zdlii{z). (5.13) 

-oo Jo 

Using the Rodrigues formula (2.18) for Pn,m on the right hand side we get 

w{z) [Pn,m{z)fi{z) - qn,m{z)] dz = C z^ log 2;D" i [z^ D'^-^ {{qz] q)n+mz"')] dfii{z), 



with C given by expression (4.6). Repeated application of summation by parts (2.11), 
gives 



W{z)[pn,miz)fi{z) - qn,miz)] dz 

-oo 

= i-irq'^C £ D^^iz"' log z]z"D^-, [{qz;qUmZn d^i,{z). 
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For D"; [z"^ log 2;] we use the g-analogue of Leibniz' rule (2.12) to find 



D-^[z"^\ogz] = Y: 



k=0 



K-'-"')l„,..,K'°. 



:x 



x=q''Z 



(5.14) 



Equation (4.7) gives us an expression for D^ z™ when n — k < m. For n — k > m we 
know that D""^^;"^ is 0. It is also easy to see that for A; > 



D'logx = {~l)'q 



k -k(k-l)/ 



l)/2 log 9 {<i\<l)k-l 



[1 — qY X 



k ' 



so it turns out that 



W{z)\pn,rn{z)h{z) - qn,m{z)\ dz = (-1) V^ ^ 



iq;q)m (ij 



m—n+k 



k=n—m L J g 
1 



,{l-qY{q]q)m-n 



+k 



X {-lfq-^^^-^)l\q-q),_^\ogqj^ z^D^^-, [(g^; g)„+„z™] dfi^{z). (5.15) 

We now have to evaluate a g-integral and a finite sum. 
(1) We start with the integral in (5.15). Apply summation by parts m times to find 

"\™D™i [{qz;qY+^z"'] dfi.iz) = (-l)™g™ |^' D™[^™](g^; g)„+„z™ rf/ii(z). 

When we use equation (4.7) for D^[z"^] and rewrite the g-integral as a infinite sum 
using (2.1), it turns out that 



z"^D"l, [(g^;g)„+„^™] dfi,{z) = (-l)™g^ 



(?;?)r 



,_,)„§(.') (.-;.). 



m+1 



We can compute the sum on the right hand side by using the g-binomial series (2.21). 
This gives 

oo • oo . f n+m+1. „\ 



3=0 



j=0 



(?; q)i 



(g; q)r 



(q 



n+2m+2. 



,q)oc 



(g-+i;g). 



{q;q)r 



(g"+-+i;g)„+i' 



So we find that 







z^D^-, [{qz;qY^^z^] df,,{z) = (-l)-g^ 



{q;q)m iq;q)m 

(l-g)™(g"+-+i;g)^+i' 



(5.16) 
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(2) Now we will work out the finite sum 



E 

k=n—m 



{q;q)m [^'') 



m—n+k 



-'g 



(1 - g)" (g; q)m-n+k 



i-l)'q-'^'-'^^\q;q)k-i 



in (5.15). Changing the index of summation gives 



E 
fe=0 



n 
k + n — m 



-'1 



(1-g)" {q;q)k 



k+n—m^—{k+n—m)(k+n—m—l)/2(„. ^\ 

q {q,q)k+n-m-i- 



Writing out the g-binomial coefficient and rewriting the powers of q gives that this 
is equal to 

" ,fc+„-^(g; g)«(g; g)r»(g; q)k+n-m-l gfc(fc+l)/2g>»«g-n(n-l)/2g-m(n^+l)/2 



E(-i) 

fc=0 



Use (g; q) 



{q;q)k+n-m{q;q)m-k{q;q)k (1 - g)" 

i/iq'i q)k+n-m = 1/(1 ~ qk+n-m^ ^^ rcwritc this sum as 



fc+n— m— 1/ y-li q Jk+n—m 



1 \n+m rnn —n{n—l) /2 —m{m+l) /2 



jq'' q)n v^/iu 
:i-a)"e.^ ' 



k=0 



m 
k 



q 



fc(fe+l)/2 



Q Qk+n-m\ ' 



Now use 



to find 



1 



1 gk+n- 



E (?'^"" 



^=0 



n+m mn —n 



I -| \n+m mn 



(n-l)/2 -m(m+l)/2_W_97n_ V^ ( n-m\^ 



v'j' 'iM y^ I n-m\ ST^ 
1 _ \n Z^ V^ y Z^ 

-^ y.' ^=0 A:=0 



,fc\^+^ fe(fe-l)/2 



g g 



If we now use the g-analog of Newton's binomial formula (2.19) with x = g^+^ then 
the expression reduces to 



/ -lyi+mmn —n{n ^)i"a 



l)/2 -m(m+l)/2_(9j_97n_ V^ { n'^-'^V (J+^- n") 



Using the g-binomial series (2.21), we have 

oo „ oo „ ( „m+\. „\ 

= (?;?)»« 

= [q] q)n-m-l 



(g"+^g)oo 

(9;g)m 



{q;q)n' 
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With this we then find 



E 



k=n—m L J g 



{q;q)m M 



m—n+k 



[I - qY {q; q)m-n+k 



{-lfq~'^'~'^'\q-,q)k-i 



= (-l)-+-g-"g-"("-l)/2g— (™+l)/2^iiil^(g;g)„_„_,. (5.17) 

Using this result together with (5.16) and the expression (4.6) for C, we find that (5.15) 
gives the required result. □ 

The integral in (5.12) can also be evaluated for m = n but with much more effort, which 
is the main reason why we have chosen m = n — 1 for our rational approximants. Lemma 
5.1 and (5.11) now imply that for m = n — 1 



\bnC,i2)-an\<di^_,ip)(q'--'] q 



-A" n2 (9;9)Li 



{<l^'';(l)n{q;q)n 



(5.1^ 



Now we can prove Theorem 1.2 

Proof of Theorem 1.2. From Section 3 we know that a„ and 6„ are integers. From (5.10) 
it follows that 6„C(?(2) — a„ 7^ because the integrand has a constant sign. Furthermore 
Lemma 3.1 gives that 



}^{dL-iiP) 



Ij-n? 



v 



24/7r2 



From(5.18) we see that 



il/n^ 24 2n-l 24 „ _3(7r^-8) 

lim o„(;o(2) — a„ < n^ lim q ^ q = p^q = p ^ 

n^oo -J ^ ' n— >oo 



< 1 



SO that lim„^oo |^nC(?(2) — a„| ^ for n -^ 00. The irrationality now follows from Lemma 
1.1. D 



5.4 Measure of irrationality for Cq{2) 



Theorem 1.2 gives rational approximations an/bn for Cq{2) that satisfy 



C,(2) 



O 



V^" 



-«(-^-«)+e)n2\ 



V 



(5.19) 



/ 



for every e > 0, with 6„ = d2n-i{p)Pn,n-i{p'^"' ^)- Combining the result of Lemma 3.1 with 
the asymptotic behavior oi Pn,n-i{p'^"'~^) given in (4.13) gives 



lim \br 



|l/n2 



lim 

n— »oo 



d2n-liP)Pn,n-l{p ) 



l/n^ 



6(7r^+4) 



p24/n^p6 ^ p-^ 
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Together with (5.19) this gives 



C,(2) - ^ 



O 






(5.20) 



, 2(,r-:+4) 



for every e > 0, which imphes the following upper bound for the measure of irrationality 

2(7r2 + 4) 37r2 



ricm < 1 



15.8369. 



TT^ 



TT^ 



6 Linear independence of 1, ^g(l), ^^(2) over 



In this section we prove Theorem 1.3 by using Lemma 1.2. To apply this lemma we 
need rational approximations for Cg(l) ^-^^d C(?(2) with common denominator. From the 
equations (4.3)-(4.4) and (5.2)-(5.3), it follows that we can take Pn = hn = d2n-i{p)Pn, 



ql = d2n-i{p)an and r* = a„, i.e. 



K = ^2n-l(pK,n-l(P " ) 



rn = 4n-liP) 



2n-2 



gn,n-l(/"-')+Pn,„-l(/"-') E 



r„,„_i(p2"-^)+Pn,„-i(p'"-^) E 



k=l P - 
2n-2 ^ 



^ ^ +(2n-l)g„,„_i(p2«-i) 



\p^ -I 



(6.1) 
(6.2) 

. (6.3) 



Then Theorems 1.1 and 1.2 imply that p* , g* and r* are integers. The following verification 
of the 3 conditions of Lemma 1.2 will establish the linear independence of 1, C(?(l) and 
Cg(2) over Q. 

6.1 Condition 1 



We will now show that ap* + 6g* + cr^ ^ for all n G N for which 2n — 1 , is prime and 
2n — 1 > c. It is sufficient to prove that 

«K + ^<i*n + cr* ^ mod $2n-i(p) 



or 



$2„-i(p) tap; + &?: + < 

for these values of n. We prove this in three steps. 



Step 1: In the first step we prove that 



«Pn + bq*n + cr*^ 



(p2n-l _ 1)2 



mod <l>2„-i(p). 



(6.4) 
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First note that d2n-i{p) = $2n-i(p)'^2n-2(p), which imphes that ^2n-i{p) divides the 
integers rfL-i(p)Pn,n-i(p^""^) and c?L-i(p)?n,n-i(p^""^) in Z. Using this we get 



aPn + bQn + crl = crf2^_i(p)r„,„_i(p2" ^) mod ^2n-\{v) 

= cd2n-l{p)d2n~2{p)^2n~l{p)rn,n^l{p'^"'~^) Hiod ^2n-l{p)- 

Expression (3.18) imphes that all the terms in the sum C(i2n-i(p)c^2n-2(p)^n,n-i(p^"~^) are 
integers except for the one with r = i = 2n — 1. So we see 



aP*n + bQn + crl = -c 



3n-2 
n — 1 



3n-l 
n 



{p2n-l_iy 



mod $2n-i(p)- 



Now we only have to eliminate the binomial numbers. This can be done by the following 
lemma. By applying it twice, with (n, m) replaced by (n, n — 1) and (n — 1, n), we obtain 
equation (6.4). 



Lemma 6.1 The following congruence for polynomials in 7j[x\ holds 

= 1 mod ^n+m{x) 



n + 2m 
m 



(6.5) 



for all n,m eN. 



Proof. We prove this result by induction on m. Obviously relation (6.5) is satisfied for 
all n G N when m = 0. Suppose that 



n + 2m — 2 
m — 1 



= 1 mod $„+m_i(a;) 



(6.6) 



for all n G N, then we can prove (6.5) as follows. The g- version of Pascal's triangle identity 
(3.16) gives that 



n + 2m 
m 



n + 2m — 1 
m 



+ X 



n+m 



n + 2m — 1 

772—1 



Note that (3.21) implies that x""*"™ = 1 mod $„+m(a;). We can also write, using (3.21), 



n + 2m — 1 
m 



yX, Xj,fi^2m—l 
[X^ X)fn\X^ X jn+m—1 



n+2m— 1 

n i^-^' 



v=n+m 



na 



X 



v=l 



n+2m—l 

n n"^^(^) 

u=n+m d\u 
m 

u=l d\u 



Since 



n+2m— 1 
m 



is a polynomial in x with integer coefficients, the cyclotomic polynomials 



$rf(a;) are irreducible over Q, and ^n+mix) is a factor of the numerator and not of the 
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denominator, it follows that ^n+m{x) divides 



n+2m— 1 
m 



n + 2m 
m 



n + 2m — 1 
m — 1 



. Using this, it turns out that 
mod ^n+m{x). 



Now we can apply the induction hypothesis (6.6) with n + 1 instead of n. This proves the 
lemma. □ 



Step 2: In the second step we prove that 

For this we need the following result of Legendre (see, e.g. [7]) 



(6.7) 



Lemma 6.2 For all positive integers p, s and every cylotomic polynomial $„ we have that 
if s I $n(p) then s = ni + 1 for some i G Nq, or s \ n. 

Suppose gcd (c, ^2n~i{p)) = s > 1. Then s \ $2n-i(p) so from Lemma 6.2 it follows that 
s = £{2n — 1) + 1 for some £ G Nq or s | 2n — 1. Since s > 1 and 2n — 1 is a prime number, 
we have in both cases that s > 2n — 1. But s \ c and we only consider values of n for 
which 2n — 1 > c. This gives a contradiction and proves equation (6.7). 



Step 3: In the last step we prove that 



*2n-l(p)t 



dln^lip) 



{P 



,2n~l 



-1 



|2- 



(6.^ 



We will do this by contraposition: suppose that there exist an integer A ^ such that 

dln-liP) 



(p2n-l _ ly 



A^2n^l{p) 



or, when we use (3.19)-(3.21), 



2n-l 



n <l>l{p)=A<^2n-lip). 



k=l 
fc|2n-l 

Suppose s is a prime number so that s \ ^2n-i{p), then it follows from the previous 
equation that there exists an integer k, with 2 < k < 2n — 2, such that s \ ^k{p)- We 
first prove that this implies gcd{^2n-i{x),^kix)) = 1 in (Z/sZ)[x]. For this we need the 
following lemma [17, Section 4.9]: 

Lemma 6.3 If s is a prime not dividing n, then x" — 1 has no repeated factors in 

{Z/sZ)[x]. 

Since s \ $2n-i(p) we can first apply Lemma 6.2 and distinguish between two cases. 
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case 1: s = i{2n — 1) + 1 for some i G Nq. In this case it is obvious that s \ 2n — 1, 
and also s \ k because A; < 2n — 2 < s. So we can apply Lemma 6.3 with parameter 
k{2n — 1) instead of n. We get that x''^'^"'~^^ — 1 has no repeated factors in (Z/sZ)[x]. 
This means, using (3.21), that 

^k{x) ■ '^2n-i{x) ■ n "^dix) 

d\k{2n-l) 

has no repeated factors, so gcd($2n-i(a^)) ^ki^)) = 1 in (Z/sZ)[x]. 

case 2: s | 2n — 1. Since s and 2n — 1 are both prime, we see that s = 2n — 1 in this case. 
Since s is prime, we also know that $s(a;) = {x'^ — l)/{x — 1). In {Z/sZ)[x] this implies 

$.(x) = ^ = ^^^ = (x-ir\ 

X — 1 X — 1 

So gcd($fc(a;), ^2n-i{x)) = gcd(<I>fe(a^)! {x — 1)'^^"'^) = (x — 1)", for some u with < m < 
s — 1 in (Z/sZ)[a;]. We know that s is a prime not dividing k, so applying Lemma 6.3 
gives that x^ — 1 has no repeated factors in (Z/sZ) [x\ . Now we can write x'^ — 1 as 

x^ - I = ^i{x) ■ ^k{x) ■ n ^d{x). 

d\k 
l<d<k 

and since <l>i(x) = x — 1, it follows that x — 1 can not be a factor of $a;(x), so we can 
conclude that gcd($2n-i(a^), $fc(a^)) = 1 in (Z/sZ)[x]. 

So in both cases we have gcd{^2n-i{x),^k{x)) = 1 in (Z/sZ)[x]. Integer division gives 
us two polynomials u{x) and v{x) such that u{x)^2n-i{x) + v{x)^k{x) = 1 in (Z/sZ)[x]. 
Evaluating this expression at x = p gives in Z the equality 

M(p)$2n-i(p) + f (p)$fc(p) = 1 + sw, 

where w is an integer. But from s \ $„+m(p) and s \ ^k{p) "we see that s | 1, which gives 
a contradiction and proves (6.8). 

From these three steps it follows that $2n-i(p) t '^Pn + ^1n + '^^n ^^r infinity many values 
of n, so the first condition of Lemma 1.2 is satisfied. 

Remark: For the case c = 0, 6 7^ we can work analogously by considering 

«Pn + bq^ + cr*n mod d2n-i{p)^2n~i{p)- 
The case c = 6 = is obvious. 
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6.2 Condition 2 and 3 

From the definition of p* and g*, it follows that |p*Cg(l) — Qnl = d2n-i{p)\PnCq{^) '~ c^nl- 
By using Theorem 1.1 together with Lemma 3.1, it turns out that 

b:c,(i) - <r/"' = (rf2„-i(p))'/"' i/3nC,(i) - «nr/"' = p-"^ < i, 

so that \PnCq{^) ~ d ~^ when n — » oo. 

It also follows from the definition of p* and r* that |p*Cg(2) — r*| = |6„^q(2) — a„| and 
from Theorem 1.2 we know already that this expression tends to zero when n tends to 
infinity. 
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